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, Abstract 

, Using the foundations laid down in Hardy and Harris [17], we present new spine proofs 

' of the /^''-convergence {p > 1) of some key 'additive' martingales for three distinct models 

_^ ' of branching diffusions, including new results for a multi-type branching Brownian motion 

. and discussion of left-most particle speeds. The spine techniques we develop give clear and 

' simple arguments in the spirit of the conceptual spine proofs found in Kyprianou [2n] and 

^ I Lyons et al [31, 30, 27], and they should also extend to more general classes of branching 

I 1 1 diffusions. Importantly, the techniques in this paper also pave the way for the path large- 

deviation results for branching diffusions found in Hardy and Harris [16, 15]. 

> : 

^ ; 1 Overview 

O ■ 

I ' In this article we use a change of measure together with spine techniques to analyze the C^- 

T-H , convergence properties (for p > 1) of the strictly-positive 'additive' martingales for three different 

■ models of branching diffusions. It is a common feature of these diffusion models, where there 
I is actually a family of such martingales {Z\ : A S M}, that for all A within an open interval 

■ about the origin the martingale Z\ is convergent in for some p > 1 subject to a suitable 
I p*"^— moment {p > 1) or LlogL {p = 1) condition on the offspring distribution; for A outside 
' of this interval, or if the L log L condition fails for the offspring distribution, the limit of Z\ is 
\ almost surely null. 

^ ■ The first model we consider is a branching Brownian motion (BBM) with random family 

sizes. After introducing the fundamental 'additive' Zx martingales and describing a 'spine' 
construction for the BBM under a change of measure using Z\, we will recall Kyprianou's [28] 
5^ I £^ -convergence result before stating necessary and sufficient conditions for ^C^-convergence of 

the Z\ martingales. Our new proof of martingale /^^-convergence for BBM uses 'spine' change of 
measure techniques and, early on in Section 2, we will include a summary of the underlying space 
and filtrations that we shall use for our spine techniques throughout this paper; a foundation 
article [17] contains full details, but we have tried to keep this article reasonably self-contained. 
Note that for BBM, the spine construction first appeared in Chauvin and Rouault [5], whilst 
Kyprianou [28] really exploited 'spine' methods in his proofs, however our spine approach does 
possess some significant differences from others. Neveu[32] used classical techniques for £P- 
convergence in the special case of binary branching. Also see Harris [21] for further discussion 
of martingale convergence in BBM and applications. 
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In Section 3, we look at a finite- type BBM model where the type of each particle controls 
the rate of fissions, the offspring distribution and the spatial diffusion. First, we will extend 
Kyprianou's [28] approach to give the analogous /^^-convergence result for this multi-type BBM 
model. We will also briefly discuss the rate of convergence of the martingales to zero and the 
speed of the spatially left-most particle within the process. Next, we give a new result on C^- 
convergence criteria, extending our earlier spine based proof developed for the single-type BBM 
case. 

The third model of Section 4 has a contrnwoMs-type-space where the type of each particle 
moves independently as an Orstein-Uhlenbeck process on R. This branching diffusion was first 
introduced in Harris and Williams ['2n] and has also been investigated in Harris [22], Git et al 
[14] and Kyprianou and Englander [7]. 

Proofs for each of these models each run along similar lines and the techniques are quite 
general, and it is a powerful feature of the spine approach that this is possible. More classical 
techniques based on the expectation semigroup are simply not able to generalize easily, since 
they often require either some a priori bounds on the semigroup or involve difRcult estimates 
- for example, in Harris and Williams [20] their important bound of a non-linear term is made 
possible only by the existence of a good theory for their operator, and this is not generally 
available. 

Briefly, to prove that the martingale converges in for some p > 1 we use Doob's theorem, 
and therefore need only to show that the martingale is bounded in C^. The spine decomposition 
is an excellent tool here for showing boundedness of the martingale since it reduces difRcult 
calculations over the whole collection of branching particles to just the single spine process. We 
find the same conditions are also necessary for ^C^-boundedness of the martingale when p > 1 by 
just considering the contributions along the spine at times of fission and observing when these 
are unbounded. Otherwise, to determine whether the martingale is merely /^^-convergence or 
has an almost-surely zero limit, we determine whether the martingale is almost-surely bounded 
or not under its own change of measure - this was Kyprianou's [28] approach and relies on a 
measure-theoretic result and has become standard in the spine methodology since the important 
work of Lyons et al [31, 30, 27]. Spine and size-biasing techniques have already proved extremely 
useful in many other branching process situations, for example, also see Athreya [2], Biggins and 
Kyprianou [3], Geiger [8, 9, 10, 11, 12], Georgii and Baake [1.3], Iksanov [24], Olofsson [33], 
Rouault and Liu [29] and Waymire and Williams [So], to name just a few. 

There are a number of reasons why we may be interested in knowing about the convergence 
of a martingale: in Neveu's original article [32] it was a means to proving /^-'^-convergence of 
martingales which can then be used to represent (non-trivial) travelling-wave solutions to the 
FKPP reaction-diffusion equation as well as in understanding the growth and spread of the 
BBM, whilst Git et al [14] and Asmussen and Bering [1] have used it to deduce the almost- 
sure rate of convergence of the martingale to its limit. Of equal importance are the techniques 
that we use here: the convergence of other additive martingales can be determined with similar 
techniques, for example, see an application to a BBM with quadratic breeding potential in 
J.W.Harris and S.C.Harris [19]; similar ideas have also been used in proving a lower bound for 
a number of problems in the large-deviations theory of branching diffusions - we have used the 
spine decomposition with Doob's submartingalc inequality to get an upper-bound for the growth 
of the martingale under the new measure which then leads to a lower-bound on the probability 
that one of the diffusing particles follows an unexpected path - see Hardy and Harris [16] for 
a spine-based proof of the large deviations principle for branching Brownian motion, and see 
Hardy and Harris [i n] for a proof of a lower bound in the model that we consider in section 4. 
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2 Branching Brownian motion 



Consider a branching Brownian motion (BBM) with constant branching rate r, which is the 
branching process whereby particles diffuse independently according to a (driftlcss) Brownian 
motion and at any moment undergo fission at a rate r to produce a random number of offspring, 
1 + A, where A is an independent random variable with distribution 

P{A = i)=p,, ie{0,l,...}, 

such that m := P{A) — 'Yli'^Lo^Pi < Offspring move off from their parents point of fission, 
and continue to evolve independently as above. We suppose that the probabilities of this process 
are {P^ : x € M} so that is a measure defined on the natural filtration {J-t)t>o such that it 
is the law of the process initiated from a single particle positioned at x. 

Suppose that the configuration of this branching Brownian motion at time t is given by the 
M- valued point process Xt := {Xu{t) : u & Nt} where Nt is the set of individuals alive at time 
t. It is well known that for any A G K, 

ueNt uGNt 

where E\ := —Xcx ■= ^A^ + rm defines a positive P-martingale, so Z\{oo) :~ limt^oo Z\{t) is 
almost surely finite under each P^. 

We are going to use a change of measure together with the so-called spine decomposition to 
determine the conditions under which this martingale is £^'(P)-convergent for some p > \. 



Theorem 2.1 If we define the measure via 



dP^ 



Zxjt) 
ZxiO) 



e 



-Xx 



Zx{t), (2) 



then it follows that under the point process evolves as follows: 

• starting from position x, the original ancestor diffuses according to a Brownian motion on 
R with drift A; 

• at an accelerated rate (1 + m)r the particle undergoes fission producing 1 + A particles, 
where the distribution of A is independent of the past motion hut is size-biased.' 

Q,(i^,)= (i±iM, ze{0,l,...}. 

m + 1 ^ ' 

• with equal probability, one of these offspring particles is selected; 

• this chosen particle repeats stochastically the behaviour of the parent with the size-biased 
offspring distribution; 

• the other particles initiate, from their birth position, an independent copy of a P' branching 
Brownian motion with branching rate r and family-size distribution given by A ( which is 
without the size-biasing). 

In this construction, the individuals that are selected to have a drift of A make up a (random) line 
of descent which has come to be referred to as the spine. The phenomena of size-biasing along 
the spine is a common feature of such measure changes when random offspring distributions 
are present. This spine construction for BBM can also be seen in Chauvin and Rouault [5], 
Kyprianou [2n] and J. Harris, S.C.Harris and Kyprianou [18]. 

In particular, Kyprianou [2s] used this change of measure and other spine techniques to 
give necessary and sufficient conditions for /^^-convergence of the Z\ martingales. By a natural 
symmetry, without loss of generality we will throughout suppose that A < 0, then: 
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Theorem 2.2 Let A ;= —\/2rm so that c\ := —E\/\ attains local maximum at X. For each 
a; e M, the limit Z\{oo) :~ limt^oo Z\(t) exists -almost surely where: 

• */ < A then Zx{oo) = -almost surely; 

• if \ <E (A, 0] and P{A\og^ A) — oo then Zx{oo) = P^ -almost surely; 

• if X & (A, 0] and P{A\og^ A) < oo then Z\{t) Z\{oo) almost surely and in C^{P^). 

We will give a generalisation of this result to a multi-type BBM in the next section, proving 
it there by extending Kyprianou's proof for BBM as well as discussing the rate of convergence to 
zero of the martingales and left-most particle speeds. In fact, in many cases where the martingale 
has a non-trivial limit, the convergence will also be much stronger than merely in C^{P^), as 
indicated by the following £f -convergence result: 

Theorem 2.3 For each x gM., and for each p G (1, 2]; 

• Z\{t) Z\{oo) almost surely and in CP{P^) if pX^ < 2mr and P{Ap) < oo 

• Zx is unbounded in Cp{P''), that is, P-^(Za(oo)p) = oo if pX^ > 2mr or P{Ap) = oo. 

We shall give a spine-based proof of this £^-convcrgence theorem in Section 2.3, but also 
see Neveu [32] for sufficient conditions in the special case of binary branching at unit rate using 
more classical techniques. Iksanov [ ' '] also uses similar spine techniques in the study of the 
branching random walk. 

2.1 The underlying space and filtrations 

For a clear understanding of the spine techniques that we shall use in all our models, we need a 
more precise description of spines than the pathwise construction given in Theorem 2.1. We give 
fuller details in Hardy and Harris [17], but to make this article self-contained we now briefly lay 
out the principal elements. The reader who is familiar with the work of Lyons et al or Kyprianou 
[28] will notice significant differences in our approach via our use of the filtrations on the single 
underlying space. 

The basic ideas of our approach is quite straightforward: given the original BBM, we first 
create an extended probability measure by enriching the process through (carefully) choosing at 
random one of the particles to be the so called spine. Now, on this enriched process, changes 
of measure can easily be applied that only affect the behaviour along the path of this single 
distinguished 'spine' particle; in our case, we add a drift to the spine's motion, increase rate 
of fission along the path of the spine and size-bias the spine's offspring distribution. However, 
projecting this new enriched and changed measure back onto the original process filtration 
(that is, without any knowledge of the distinguished spine) brings the fundamental 'additive' 
martingales into play as a Radon-Nikodym derivative. The four probability measures, various 
martingales, extra filtrations and clear process constructions afforded by our setup, together 
with some other useful properties and tricks, such as the spine decomposition, provides a very 
elegant, intuitive and powerful set of techniques for analysing the process. 

All three models that we consider in this article shall be built on the same underlying space 
of sample trees with spines, and the measures will all be constructed in analogous ways. Here 
we lay out the details for the current model of branching Brownian motion, and leave it to the 
reader to bridge the details when it comes to the other models - [17] contains all the details in 
a more abstract setting that will cover every model considered in this article. 

The set of Ulam-Harris labels is to be equated with the set ft of finite sequences of strictly- 
positive integers: 

n := {0} U U (N)", 

nGN 
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where we take N = {1,2,...}. For two words u,v € ft, uv denotes the concatenated word 
(u0 = %u = u), and therefore Q. contains elements like '213' (or '0213'), which we read as 'the 
individual being the 3rd child of the 1st child of the 2nd child of the initial ancestor 0'. For two 
labels £ the notation v <u means that v is an ancestor of u, and |u| denotes the length 
of u. The set of all ancestors of u is equally given by 

\y : V < ~ \y : 3w G such that vw = w}. 

Collections of labels, ie. subsets of fi, will therefore be groups of individuals. In particular, 
a subset r C f2 will be called a Gallon- Watson tree if: 

1. e T, 

2. if u, w € ri, then uv G t implies u € t, 

3. for all u E T, there exists Au G 0,1,2,... such that uj € t if and only ii 1 < j < 1 + Au, 
(where j G N) . 

In general, each node u of the underlying trees will have l + Au branches coming from it, whereas 
the special case of binary-branching, for example, would correspond to A„ = 1 at every node. 

The set of all Galton- Watson trees will be called T. Typically we use the name r for a 
particular tree, and whenever possible we will use the letters u or w or w to refer to the labels 
in r, which we may also refer to as nodes of r or individuals in r or just as particles. 

Each individual should have a location in M at each moment of its lifetime. Since a Galton- 
Watson tree r € T in itself can express only the family structure of the individuals in our 
branching process, in order to give them these extra features we suppose that each individual 
u G T has a mark (X^, (t„) associated with it which we read as: 

• <7u G M+ is the lifetime of m, which determines the fission time of particle u as Su ■= 
Ei)<u '^v (with 5*0 := (T0). The times Su may also be referred to as the death times; 

• Xu : [Su — ctm, Su) ^ M gives the location of u at time t € [Su — ctu, Su). 

To avoid ambiguity, it is always necessary to decide whether a particle is in existence or not at 
its death time. 

Remark 2.4 Our convention throughout will be that a particle u dies 'just before' its death time 
Su (which explains why we have defined Xu : [Su — cTu^Su) ■ for example). Thus at the time 
Su the particle u has disappeared, replaced by its 1 + Au children which are all alive and ready 
to go. 

We denote a single marked tree by (r, X, a) or (r, M) for shorthand, and the set of all marked 
Galton- Watson trees by T: 

• T := |(r,X,(T) : r G T and for each u G r, (t„ G M+,X„ : [Su - au,Su) 

• For each (r, X, a) G T, the set of particles that are alive at time t is defined as Nt := {u G 
T : Su — <Ju <t<Su). 

For any given marked tree (r, M) G T we can identify distinguished lines of descent from the 
initial ancestor: 0,ui,U2,U3, . . . G r, in which U3 is a child of U2, which itself is a child of ui 
which is a child of the original ancestor 0. We'll call such a subset of r a spine, and will refer to 
it as ^: 
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• a spine ^ is a subset of nodes {0, ui, U2, U3, . . .} in the tree r that make up a unique line 
of descent. We use to refer to the unique node in ^ that that is alive at time t. 

In a more formal definition, which can for example be found in the paper by Rouault and Liu 
[29], a spine is thought of as a point on dr the boundary of the tree - in fact the boundary 
is defined as the set of all infinite lines of descent. This explains the notation ^ e 9r in the 
following definition: we augment the space T of marked trees to become 

• f {(r,A/,0 : (r, M) G T and ^ G 5t} is the set of marked trees with distinguished 
spines. 

It is natural to speak of the position of the spine at time t which think of just as the position of 
the unique node that is in the spine and alive at time t: 

• we define the time-i position of the spine as Xu(t), where u E ^ O Nt- 

By using the notation to refer to both the node in the tree and that node's spatial position 
we are introducing potential ambiguity, but in practice the context will make clear which we 
intend. However, in case of needing to emphasize, we shall give the node a longer name: 

• nodet((r, M, ^)) := u if u G ^ is the node in the spine alive at time t, 

which may also be written as nodet(^). 

As the spine diffuses, at the fission times Su for m G ^ it gives birth to some offspring, 
one of which continues the spine whilst the others go off to create subtrees like copies of the 
BBM. These times on the spine are especially important for the later spine decomposition of 
the martingale Z\, and we therefore give them a name: 

• the sequence of random times {S'„ : u G ^} arc known as the fission times on the spine; 

Finally, it will later be important to know how many fission times there have been in the spine, 
or what is the same, to know which generation of the family tree the node S,t is in (where the 
original ancestor is considered to be the 0th generation) 

Definition 2.5 We define the counting function 

nt = \ nodet{£,)\, 

which tells us which generation the spine node is in, or equivalently how many fission times there 
have been on the spine. For example, if S,t = (0,mi,U2) then both and ui have died and so 
nt = 2. 

The collection of all marked trees with a distinguished spine is T; on this space we define four 
filtrations of key importance that encapsulate different knowledge, but see Hardy and Harris [17] 
for more precise details: 

• Tt knows everything that has happened to all the branching particles up to the time i, but 
does not know which one is the spine; 

• Tt knows everything that Tt knows and also knows which line of descent is the spine (it 
is in fact the finest filtration); 

• Qt knows only about the spine's motion in J up to time t, but does not actually know 
which line of descent in the family tree makes up the spine; 
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• Qt knows about the spine's motion and also knows which nodes it is composed of. Further- 
more it knows about the fission times of these nodes and how many children were born at 
each time. 



Having now defined the underlying space for our probabilities, we remind ourselves of the prob- 
ability measures: 

Definition 2.6 For each x € R, let be the measure on (T^Too) such that the filtered prob- 
ability space {'T,!Foo^ {^t)t>o, P ) makes the M.-valued point process Xt = {Xu{t) : u G Nt} the 
canonical model for BBM. 

Our spine approach relies first on building a measure under which the spine is a single 
genealogical line of descent chosen from the underlying tree. If we are given a sample tree (t, M) 
for the branching process, it is easy to verify that, if at each fission we make a uniform choice 
amongst the offspring to decide which line of descent continues the spine ^, when u £ t we have 

Prob(riee)= HlTX- 

This simple observation at (3) is the key to our method for extending the measures, and for this 
we make use of the following representation found in Lyons [M)]. 

Theorem 2.7 /// is a J- 1 -measurable function then we can write: 

/ = E /«1(6=«) (4) 

where /„ is J- 1- measurable. 

We use this representation to extend the measures . 

Definition 2.8 Given the measure P^ on (T, JFqo) we extend it to the probability measure P^ 
on (T^Tao) by defining 

X/di^^-XE/"^T^dP^ (5) 

for each f G mj^t with representation like (4). 

The previous approach to spines, exemplified in Lyons [30], used the idea of fibres to get a 
measure analogous to our P that could measure the spine. However, a perceived weakness in 
this approach was that the corresponding measure had a time dependent total mass and could 
not be normalized to become a probability measure with an intuitive construction, unlike our 
P. Our new idea of using the down-weighting term of (3) in the definition of P is crucial in 
ensuring that we get a probability measure, and leads to the very useful situation in which all 
measure changes in our formulation are carried out by martingales. 

Theorem 2.9 This measure P^ really is an extension of P^ in that P = P\j^aa- 

The spine diffusion S,t is J-j-measurable, and it is immediate that, under P^ , the spine diffusion 
is a Brownian motion that starts at a;. In fact, it is easy to see that 

Theorem 2.10 Under P^ , the process "Kt evolves as follows: 

• starting from x, the spine diffuses according to a (driftless) Brownian motion on M,' 
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• at rate r the spine undergoes fission producing 1 + A particles, where A is independent of 
the spine's motion with distribution {pk '■ k>0}; 

• with equal probability, one of the spine 's offspring particles is selected to continue the path 
of the spine, repeating stochastically the behaviour of its parent; 

• the other particles initiate, from their birth position, independent copies of a P branching 
Brownian motion with branching rate r and family-size distribution also given by A. 

It is useful to have this natural construction of P in mind when we start to change measure. 
2.2 New measures for BBM 

Having seen the construction of the underlying space and the measure P^, we can define a 
measure Qa via a Radon-Nikodym derivative with respect to P after recalling three simpler 
changes of measures. 

Changing measure with the P-martingale 



would make the spine process a Brownian motion with drift A under the new measure. 
We also note that 

e-""^*(l + m)"' 

is a P-martingale that will increase the rate of the Poisson process Ut of fission times on the 
spine from r to (1 + m)r after a change of measure. 
Lastly, 

^\ \ l + m 

is also a P-martingale that will produce size-biasing (only) along the spine, giving an offspring 
distribution {(1 + k)pk/{l + m) ■.k>0}. 

Combining these components into a single change of measure leads to the following key result: 

Theorem 2.11 Define the measure Q\ on (T, J^oo) by: 



dP^ 



Under Q^, the process evolves as follows: 

• starting from x, the spine ^4 diffuses according to a Brownian motion with drift A on R; 

• at accelerated rate (l-\-m)r the spine undergoes fission producing l-\-A particles, where A is 
independent of the spine's motion with size-biased distribution {{l-\-k)pk/ {l-\-m) : k > 0}; 

• with equal probability, one of the spine 's offspring particles is selected to continue the path 
of the spine, repeating stochastically the behaviour of its parent; 

• the other particles initiate, from their birth position, independent copies of a P branching 
Brownian motion with branching rate r and family-size distribution given by A, that is, 
{pk : fc > 0}. 

The measure Qx that we introduced in Theorem 2.1 via its pathwise construction can equiv- 
alently be obtained from Q\ by restricting it to the filtration {J^t)t>i3 on which the original 
measure P is defined, that is, simply ignoring information identifying the spine: 
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Theorem 2.12 If we define := Q^l^^oo, then < 



IS a measure on 



Tea with 



(7) 



Under Q^, the branching- diffusion point process can be pathwise constructed exactly as de- 
scribed in Theorem 2.1. 

There are at least two ways to prove this resuh: Kyprianou [28] bases his proof on a decomposi- 
tion of a related (non-probability) measure P* as a product of measures for the spine's motion, 
the fission-counting process nt, and measures on the sub-trees born from the spine. Since we 
have clear constrictions of the probability measures P and Qa and richer collection of filtrations. 
we have an alternative: 

Proof of Theorem 2.12: Since <Ql := Q^|^<^, whilst ^ P'-']^^ from Theorem 2.9 and 
it is clear that the change of measure at (6) projects onto the sub-algebra as a conditional 
expectation, we have 



dP^ 



dP^ 



-Xx 



Using representation (4) and recalhng from (3) that P{^t = n\Tt) = OtKu (1 + ^u) ^ yields 



E 



[](i + A„)xP(6 = u|.Ft; 



^XX^{t)-E,t ^ z^{t) 

u<£Nt 



and trivially noting Z\{0) ~ e under P^ completes the proof. 



□ 



2.3 Proof of Theorem 2.3 

Just before we proceed to the proof we recall the naturally occurring eigenvalue E\ := ^A^ +mr, 
noting that under the symmetry assumption that A < and for p G (1,2]: 

pEx - Epx > ^ cx> Cpx <^ pX^ < 2mr 

and that this always holds for some p > 1 whenever A g (A, 0], that is, when A lies between the 
minimum of cx found at A and the origin. 

2.3.1 Proof of part 1: 

We are going to prove that for every p G (1,2] the martingale Zx is £P(P)-convergent if pEx — 
Epx > 0. Furthermore, since P^{Zx{t)P) — e^^^ P^ [Z x{tY) we do not lose generality supposing 
that a; = 0; from now on this is implicit if we drop the superscript by simply writing P. 
From the change of measure in Theorem 2.1 or Theorem 2.12 it is clear that 

p{Zxm = Pizxity-'zxit)) = ^ixizxm. 

where q := p — 1. Our aim is to prove that Q^xiZxiiY) is bounded in i, since then Z^{t) must 
be bounded in ^{P) and Doob's theorem will then imply that Zx is convergent in Cp{P). 
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As we mentioned, the algebra Qoo gives us the very important spine- decomposition of the 
martingale Zx: 

Qx{Zx{t)\goo) = ^ +e^«*-^-*, (8) 

fc=i 

where Ak is the number of new particles produced from the fission at time Sk along the path 
of the spine, and the sum is taken to equal if = 0. Full details of this can be found in 
[17], but the intuition is quite clear: since the particles that do not make up the spine grow to 
become independent copies of distributed as if under P, the fact that Zx is a P-martingalc 
on these subtrees implies that their contributions to the above decomposition are just equal to 
their immediate contribution on being born at time Sk at location ^5^. . Remark: We emphasize 
that here we must use the measure Qa, since Qa cannot measure the algebra Goo ^ J^oo- 

We can now use the conditional form of Jensen's inequality followed by the spine decompo- 
sition of (8) coupled with the simple inequality. 

Proposition 2.13 If q e (0, 1] and u,v > then {u + w)' < u'^ + v'^ , 
to obtain, 

^x{Zxmgoo)<Qx{Zx{t)\gooy (9) 

< ^|e«^««'=-«-^^'^'' + e«^«'~«^^*. (10) 
fc=i 

With the tower property of conditional expectations and noting that Qa and Qa agree on Tt, 
^x{Zx{tY) = ^x{Zxm - Qa [^x{ZxmQoo)) (11) 

<Qa($3^^'^^^-~«''^^'=) +QA(e«^«'"''^"*), (12) 

k=l 

and the proof of £^'(P)-boundcdness will be complete once we show this is bounded in t. 

As written, (12) is made up of two terms, and since they play a central role from here on 
we name them explicitly: on the far right we have the spine term Qa (e'^'^^'^'''^^*) , the other 
being the sum term QA(Efcii A^e'^^^, -<?i5ASfc j ^ 

The spine term: Changing from P to Qa gives the spine a drift of A, and therefore the change- 
of-measure for just the spine's motion (i.e. on the algebra Qt) is carried out by the martingale 
gAjt-jA^t^ so 

since the second-line term e^''^^*"^'^^^"* is also a P-martingalc and j{p\Y — ^A^ = Epx — Ex- 

The sum term: Recall, under the measure Qa we know that the fission times {Sk : fc > 0} on 
the spine occur as a Poisson process of rate (1-1- m)r with the fc*^ fission yielding an additional 
Ak offspring, each Ak being an independent copy of A which has the size-biased distribution 

{{l + k)pk/{l + m):k>Q}. 
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First, conditioning on the motion of the spine (without knowledge of the fission times or 
family sizes) and appealing to intuitive results from Poisson process theory (see [25] for example) 
yields 

"t . ft 

ix\ 

fe=i 

Taking expectations of both sides of (14) and using Fubini's theorem then gives 

QA(XlA^e«^«^^'=-'^^'^'=) = (l + m)rQA(^') / Qx^e"^^"-"^^^'] ds 
fe=i 



i — 1 -^0 



(l + m)rQA(i') / e-'^P^^-^">-^' ds, using (13). 







Thus we have found an explicit upper-bound (if pE\ ^ Ep\): 



px 



(15) 



\pEx - Epx L J ^ ^ 

Finally, we also observe that 

Lemma 2.14 If p G (1,2] and q:=p-l, Qa(^') < oo tf and only if P{Ap) < oo 
since 

^ — ' TO + 1 TO + 1 V'^ — ' ^ — ' / TO + 1 TO + 1 

i—1 i—1 2—1 

Hence, if we have pEx — Epx > in addition to P(Ap) < oo, this imphcs that P^'{Zx{t)P) 
will remain bounded as t — > cx), which together with Doob's theorem will complete the proof of 
the first part of Theorem 2.3. □ 



2.3.2 Proof of Part 2: 

We seek to show that Zx is unbounded in Cp{P''^) if cither pEx — Epx < or P{Ap) ~ oo. 
Note that if Zx is £^'(P^) bounded then 

P^(Za(oo)p) = hm P'-'iZxm < oo 

t — >OG 

hence, Qa(2'a(oo)'?) < oo and Za(oo)'^ is a uniformly integrable Q^-submartingale. In particular, 
for any stopping time T, (Za(oo)9|.Ft) > ZxiT)" hence Q^(Za(oo)«) > Q^(Za(T)9). 

First, by considering only the contribution of the spine Zx{t) > e'*'^*"^^* for all t > and 
recalling (13), we see that 

Ql{Zx{ty) > Q;j(e*^«*"'?^^*) = e'^^^-ipBx-E^^)t 

and Zx is therefore unbounded in Cp{P^) if pEx — Epx < 0. 

Now, let T be any fission time along the path of the spine, then 

Zx{T) > (l + i)e^«^-^^'^ 

where A is the number of additional offspring produced at the time of fission. Then, 

QliZxiT)") > Ql ((1 + e"^- Q^ACe-^^^^-^"^^^) 
and so Zx is unbounded in /:p(P^) if Ql ((1 + A)A = oo P{Ap) = oo. □ 
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3 A typed branching diffusion 



Wc move on to consider a typed branching diffusion where the type of each particle evolves 
as a Markov chain and influences the fission rate, offspring distribution and spatial diffusion 
coefficient of the particle. We will follow a similar notation and setup as for the BBM ease, but 
leave some details to the reader. 

A 'typical' particle motion will be like the process (Xt, lt)t>o in J R x / where: 

(i) the type location, Yt, is an irreducible, time-reversible Markov chain on the finite type-space 
/ := {!,..., 71} with Q-matrix 9Q (9 is a strictly positive constant) and invariant measure 
TT = (tti, . . . ,7r„); 

(ii) the spatial location, Xt, moves as a driftless Brownian motion on R with diffusion coefficient 
a{y) > whenever Yt is in state y, that is, 

dXt — a{Yt)^ dBt, where Bt a Brownian motion. (16) 
The formal generator of this process {Xt,Yt) is therefore: 

nF{x,y)^^a{y)^+eY,Qiy,j)Fi^:j), (F-.J^R). (17) 

We often use matrix calculations, and it is convenient to gather the diffusion coefficients together 
in a diagonal matrix A := diag[a(l), . . . , a{n)]. 

The typed branching diffusion. Consider a typed branching Brownian motion whereby 
individual particles move independently according to the 'typical' particle motion, as described 
above, and at any moment a particle currently of type y will undergo fission at rate r(j/) to be 
replaced by a random number of offspring, 1 + A{y), where A{y) £ {0, 1,2,.. .} is an independent 
RV with distribution 

P{A{y)=i) ^p,{y), ie{0,l,...}, 

and mean m{y) := P(^A{y)^ < 00 for all y & I. At birth, offspring inherit the parent's spatial and 
type positions and then move off independently, repeating stochastically the parent's behaviour, 
and so on. As before, we let N{t) be the set of particles alive at time t using the Ulam-Harris 
labelling, where a particle u S N{t) has spatial position and type given by {Xu{t),Yu{t)). We 
gather together the birth rates in a diagonal matrix R := diag[r(l), . . . ,r(n)] and the mean 
number of offspring in M := diag[m(l), . . . , m{n)]. 

We suppose that the probabilities for this process are given by {p^^^ : {x,y) € j} defined 
on the natural filtration {!Ft)t>o, where P^ "^ is the law of the process starting with one initial 
particle of type y at spatial position x. 

The typed branching diffusion with spine. As seen before, we can extend the measures 
{P'^ y : {x,y) £ J} to {P'^ y : {x,y) € j] by identifying a single distinguished infinite line 
of descent starting from the initial single particle, known as the spine. The natural filtration 
{^t)t>o then contains all information about the process, including the identity of the spine. 

Theorem 3.1 Under P^^y^ the process Xt := {{Xu{t),Yu{t)) : t > 0} evolves as follows: 

• starting from (x, y), the spine (^t, rjt) evolves as a process with type rjt moving as a Markov 
chain on I with Q-matrix 6Q and spatial position ^t diffusing as a (driftless) Brownian on 
R with variance coefficient a{rit); 

• if the spine is currently of type y, at rate r{y) the spine undergoes fission producing \-\-A{y) 
particles, where A{y) is an independent random variable with distribution {pkiu) ■ k > 0}; 
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• with equal probability, one of the spine 's offspring particles is selected to continue the path 
of the spine, repeating stochastically the behaviour of its parent; 

• the other particles initiate, from their birth position, independent copies of a P '' typed 
branching Brownian motion as described above. 

It should be noted that the condition of time-reversibihty on the Markov chain is not abso- 
lutely necessary, and is really just a simplifying assumption that gives us an easier theory 
for the matrices and eigenvectors; our aim is really to show how the spine techniques work ~ 
lessening the geometric complexity of the model serves a good purpose. 

Note, the special case of the 2-typc BBM model was considered in Champncys et al [4] by 
different means. Also, in our model, at the time of fission a typc-y individual can produce only 
type-y offspring. This is not the same as the case in which a type-y individual may produce 
a random collection of particles of different types - as considered in T.E. Harris's classic text 
[23], for example. Other forms of typed branching processes have also been dealt with by spine 
techniques, for example, see Lyons et al [27] or Athreya [2] for discrete-time models in which a 
particle's type does not change during its life but a type-w individual can give offspring of any 
type according to some distribution. See also the remarkable work of Georgii and Baake [13] 
that uses spine techniques to study ancestral type behaviour in a continuous time branching 
Markov chain where particles can give birth to across all types. In principle, our spine methods 
will be robust enough to extend to all these other type behaviours (with added spatial diffusion). 

3.1 The martingale 

Via the many-to-one lemma (see [17]) or generators it is easy to see that for any A e M, any 
function (vector) : / — > K and any number E\ G M, the expression 

ZA(t):= «A(rfe(t))e^^'=(*'-^^*, (18) 

«GJV(t) 

will be a martingale if and only if v\ and E\ satisfy: 

{^\^A + BQ + MRyx^Exvx, (19) 

that is, v\ must be an eigenvector of the matrix \\^A + 6Q + MR, with eigenvalue Ex. 
Definition 3.2 For two vectors u, v on I , we define 

n 

1=1 

which gives us a Hilbert space which we refer to as C^iir). We suppose that the eigenvector vx 
is normalized so that \\vx\\^ '■= {vx,vx}^ = 1- 

The fact that the Markov chain is time-reversible implies that the matrix ^X^A + 9Q + MR is 
self-adjoint with respect to this inner product. This in itself is enough to guarantee the existence 
of eigenvectors in £^(7r), but the fact that we are dealing with a finite-state Markov chain means 
that we also have the Perron- Frobenius theory to hand, which allows us to suppose that vx is 
a strictly positive eigenvector whose eigenvalue Ex is real and the farthest to the right of all 
the other eigenvalues - see Seneta [34] for details. This implies a useful representation for the 
eigenvalue: 
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Theorem 3.3 

Ex= sup {({Xy2)A + eQ + MR)v,v) , (20) 
since it is the rightmost eigenvalue. 

A proof can be found in Krcyzig [2(i]. From this it is not difficult to show that Ex is a strictly- 
convex function of A. Interestingly, it will be seen in our proofs that it is the geometry of 
the eigenvalue E\ that determines the interval that gives rise to martingales Z\{t) that arc 
/^''-convergent. 

Corollary 3.4 As a function of X, E\ is strictly- convex and infinitely differ entiable with 

E'^ = \{Avx,vx)^. (21) 

// we define the speed function 

cx := -Ex/\, (22) 

then on (— oo, 0) the function cx has just one minimum at a single point X{d) , either side of which 
cx is strictly increasing to +oo as either X J, —oo or A t 0. In particular, for each X G (A(0),O] 
there is some p > 1 such that cx > Cpx; on the other hand, if X < X{6) there is no such p > 1. 

We refer to the function cx as the speed function since it relates to the asymptotic speed of the 
travelling waves associated with the martingale Zx{t); sec Harris [21] or Champneys et al [4] for 
details of the relationship between branching-diffusion martingales and travelling waves. 

Since Zx{t) is a strictly-positive martingale it is immediate that Zx{oo) := liuit^oa Zx{t) 
exists and is finite almost-surcly under P^^y. As before, by symmetry we shall assume that 
A < and, without loss of generality, we also suppose that P{A(y) = 0) = 1 whenever r(y) = 
to simplify statements. We shall prove necessary and sufficient conditions for £^ -convergence of 
the Zx martingales: 

Theorem 3.5 For each x CzM., the limit Zx{oo) := limt^oo Zx{t) exists P^'^-a.s. where: 
' "if ^ l£ A(0) then Zx{oo) = P^'^ -almost surely; 

• tfXe {X{0),0] and P {A{y) \og+ A{y)) = oo for some y £ I, then Zx{oo) = P^-^'-a.s.; 

• if X e {X{6),0] and P {A{y) log'^ A{y)) < oo for all y e I, then Zx{t) Zx{oo) almost 
surely and in C^^P^'^). 

Once again, in many cases where the martingale has a non-trivial limit, the convergence will 
be much stronger than merely in £^(P^'^), as indicated by the following new £''-convergence 
result that we will prove by extending our earlier new spine approach: 

Theorem 3.6 For each x £ R, and for each p E (1, 2]; 

• Zx{t) Zx{oo) a.s. and in £P(P^.a) ifpEx - Epx > and P{A{y)P) < oo for all ye I. 

• Zx is unbounded in 0'{P^''^), that is, P^'^(Za(oo)^) ^ oo if either pEx — Epx < or 
P{A[y)P) = oo for some y G /. 

Note, when X < 0, the inequality pEx — Epx > is equivalent to cx > Cpx and holds for some 
p e (1, 2] if and only if X E {X{9), 0]. 



14 



3.2 New measures for the typed BBM 

With the construction of the underlying space and the measure P^'V, we can define a measure 
Qa via a Radon-Nikodym derivative with respect to P by combining three simpler changes of 
measures that only affect behaviour along the spine. 
We first alter the motion of the spine: 

Lemma 3.7 

is a P -martingale and using this martingale as a Radon-Nikodym derivative would give a new 
measure under which S,t would have instantaneous drift a{r]t)X and rjt would have modified Q- 
matrix Q\ with invariant measure ttx = v\tt . 

ft is easy to see this first martingale property with classical 'one-particle' calculations, for 
example, using the Feynman-Kac formula and noting the relation (19). We will discuss the 
behaviour of the spine under such a change of measure in detail in the next subsection. 

Secondly, we alter the fission rate along the spine: 

Lemma 3.8 

is a P -martingale that will increase the rate at which fission times occur on the spine from R{r]t) 
to {l + m{r]t))Rivt)- 

Finally, we size-bias the offspring distribution along the spine: 

Lemma 3.9 

-pr 

}X l-t-mfej 

is a P -martingale that will cause the family distribution on the spine to he size-biased to the 
distribution 

ProbiAM^^)^^l}I!^, ^e{0,l,...}. 
m(w) + 1 

Combining these components in one change of measure leads to the following key definition: 
Definition 3.10 For each A e R we define a measure Q^'^ on (T, JFqo) via 



dP^,y 



MVt) „A(6-.)-£.t -Q ^ (23) 



2A 



e 



As in the previous BBM case, a proof using the conditional expectation of this measure-change 
martingale confirms: 

Theorem 3.11 If we define Q^''^ :~ Q^''^\j^^, then Q^'^ is a measure on Too and 



dp^,y 



Notice that with this result, starting with the three simple known martingales, we have 
actually shown that Z\ must, in fact, be a martingale. This route offers a simple way of getting 
very general 'additive' martingales for the branching process. 
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3.3 The spine process {^t, Vt) under Qx 

In the BBM model it was clear to sec that the spine received a drift under the measure Qa; 
something similar happens here: 



Lemma 3.12 Under <Q)x the spine process {(.t,Vt) has generator: 

1 , .d^F , ^dF 
2a(y)^+a(y)A^ 



1 0^ F dF 
HxF [x, y) := 7^a{y)j^ + a{y)\— + ^ eQx{y.j)F{x,j), (24) 



where Qx is an honest Q-matrix: 



n ( 1 



i,i) + ^a{i) - Ex+r{i) if i ^ j 



with invariant measure ttx = w^tt. 



Thus under Qx the Q-matrix (generator) of r/t is changed to 6Qx, and the process € R is 
given an instantaneous drift of a{rit)X. The form of this above generator can be obtained from 
the theory of Doob's /i-transforms, due to the fact that on the algebra Qt the change of measure 
is given by: 

^,^(^^)e/o*M«(';.)dSgA6-s,t^ (25) 



tx 



vx{y)e 



The long-term behaviour under of the spine diffusion can now be retrieved from the 
generator (24) and the properties of Ex stated in Lemma 3.4: 

Corollary 3.13 Almost surely under Q^'^, the long-term drift of the spine is given explicitly as 

lim t-^^t = E'x 



and hence 



(.+cxt-.h (26) 

^ I -GO tfx<x{e) ^ ' 



whereas, if X ^ X the process + cxt will be recurrent on R under Qx- 
Proof: From the generator stated at (24) we can write: 



^t^Bl^J a{ris)dsj + Xj a{T],)ds, 

where B{t) is a QA-Brownian motion. Then by the ergodic theorem and the fact that ttx = x'f tt: 
^ A^a(y)7rA(2/) = A ^ a(y)z;2 (y)7r(j/) = X{Avx,vx)^ = E'^. 

yel y£l 

Direct calculation from (22) gives E'^ ~ — ca — Xc'^, and therefore t^^{£^t +cxt) ^Xc\, whence 
whether we are to the left or right of the local minimum of ca found at A determines the behaviour 
of -|- cxt, as is required. Lastly, when A = A, with the laws of the iterated logarithm in mind, 
it is not difhcult to see that both i?(/Q a{r]s) ds) and {Xa{ris) + cx) ds will fluctuate about the 
origin, hence + cxt will be recurrent under Qa- D 
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3.4 Construction of the process under Qx 

Drawing together the elements from this section, we now present the important intuitive pathwise 
construction of the new measure Q^'^: 

Theorem 3.14 Under Q^'^ , the process \t evolves as follows: 

• starting from {x,y), the spine (5t,?7t) evolves as a Markov process with generator Tlx, that 
is, rjt evolves as Markov chain on I with Q-matrix 6Q\ and moves as a Brownian motion 
on M with variance coefficient a{rit) and drift a{rjt)\. 

• whenever the type of the spine rj is in state y Cz I , the spine undergoes fission at an 
accelerated rate {l + m{y))r, producing l + A{y) particles where A{y) is independent of the 
spine's motion with size-biased distribution {(1 + k)pk{y) / {1 + nn{y)) : k > 0}; 

• with equal probability, one of the spine 's offspring particles is selected to continue the path 
of the spine, repeating stochastically the behaviour of its parent; 

• the other particles initiate, from their birth position, independent copies of P '' typed 
branching Brownian motions. 

3.5 Proof of Theorem 3.5 

The following proof is an extension of that given for BBM by Kyprianou [2n]. The second part 
of the following theorem is the key element in using the measure change (3.10) to determine 
properties of the martingale Z\: 

Theorem 3.15 Suppose that P and Q are two probability measures on a space {p,,J-oo) with 
filtration {Tt)t>o, such that for some positive martingale Zt, 



mi 
dP 



The limit Z^q '■= limsup(_,g^ Zt therefore exists and is finite almost surely under P. Further- 
more, for any F G J-qo 



and consequently 



Q{F)^ / Zoo dP + Q(Fn {Zoo (27) 



F 



(a) P{Zoo = 0) = 1 ^ Q{Zoo = oo) = 1 (28) 
(6) P(Zoo) = 1 ^ Q{Zoo < oo) = 1 (29) 

A proof of the decomposition (27) can be found in Durrett [G], at page 241. 

Suppose that A < A < 0. Ignoring all contributions except for the spine, it is immediate that 

Zx{t) - J2 «A(i;.(0)e^^"W-^^* > «A(^t)e^(«'+=^*) 

ueNt 

where, from Corollary 3.13, under the measure Q\ the spine satisfies \im.mi{^t + cxt} ~ ^oo a.s. 
and v\ > 0, hence limsup^^oo Zx(t,) = oo almost surely under Qx, yielding P(Za(oo) = 0) = 1. 

Note that, for y G I, P{A{y) log+ ^) < oo J2k>i -P(log"'" A{y) > cfc) < oo for any 

c > 0, where recall that A{y) has the size-biased distribution {{i + l)pk{y)/{l + rn{y)) : k > 0}. 
Then for an IID sequence {A„(j/)} of copies oi A{y), Borel-Cantelli reveals that, P almost surely, 

r -1, Jo if P(^(y)log+A(2/)) <cx), 

hmsupji log An{y)^< -rp/./N, + ^^ (30) 
t^oo CO if P(^(y) log^ A(?/)) = oo. 
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Now suppose that A £ (A, 0] and P{A{y) log^ ^iu)) = oo for some y & I (with r(y) > 0). Let 
Sk be the time of the k*"^ fission along the spine producing Akirjs^) additional particles, then 

where (^t + c\t)/t — Ac^ > 0, 7]t is ergodic so the event {rjs,, = y} will occur for infinitely 
many k since r{y) > 0, and rit/t Rv\,v\ >7r so Sk/k Rv\,v\ >~^, hence the super- 
exponential growth for Ak{y) from (30) gives limsupf^o^ ^^(t) = oo Q^-almost surely which 
then implies that P{Z\{oo) = 0) = 1. 

Finally, suppose that A S (A,0] and P {A{y) \og^ A{y)) < oo for all y I. Recall the very 
important spine-decomposition of the martingale Z\ from (8): 

rit 

fe=i 

In this case, the facts that {£,t+c\t)/t —t — Ac^ > and Sk/k —f< Rv\,v\ together with the 
moment conditions and (30) implying that the Ak{yys all have sub-exponential growth means 
that 

limsupQA(^A(i)|aoo) QA-a.s. 

t — ^oo 

Fatou's lemma then gives liminft_,oo ■Z^A(i) < oo, QA-a.s., hence also Q^-a.s. In addition, since 
Z\(t)~^ is a positive Q^-martingale (recall Theorem 3.11) with an almost sure limit, this means 
that limt^oo Z\{t) < oo. Q^-a.s. and then (29) yields that P(Za(oo)) = 1 and so Z\{t) converges 
almost surely and in £^(P). □ 



3.5.1 Discussion of rate of convergence to zero and left-most particle speed. 

Alternatively, when A < A we can readily obtain the rate of convergence to zero with the following 
simple argument, adapted from Git et al [14]. By Proposition 2.13, 

Zx{ty < J2 i'A(>"«(i))''e«^''^"**^'^''"^e9^'''"'''-'* < A' Z,A(i)e''^^''~''''^* 

where K := maxy^i v'^{y) / Vq\(y) < oo since / is finite and > 0. Recall that ca has a minimum 
over A € (— oo,0] at A with cj^ = = < ^''^A'^a Then, since Zq\{t) is a convergent 
martingale, we can choose q such that q\ = \ giving Z\{t) decaying exponentially to zero at 
least at rate A(ca — Cx)- 

Further, once we know that P and Qa are equivalent for every A € (A,0], since the spine 
moves such that ^j/t —c\ — Xc'-^ under Qa, the left-most particle L{t) := inf„g7v(t) ^ti(i) 
must satisfy liminff L(t)/t < — cj^, P-a.s. On the other hand, the convergence of the Z\ P- 
martingales quickly gives the same upper bound on the fastest speed of any particle, leading to 
L{t)/t —c^i P-a.s. This result also reveals that the rate of exponential decay found above is 
actually best possible. 

3.6 Proof of Theorem 3.6 
3.6.1 Proof of Part 1: 

Suppose p £ (1, 2], then with g := p — 1 a slight modification of the BBM proof arrives at 
P-'^(ZA(t)f) = e^-vxiyWx^Zxm 

< e^^vxiy) (^^^''[f^^Alvxivs.re^^i^^-''''^'^^^ 
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and the proof of £P-boundedness will be complete once we show that this RHS expectation is 
bounded in t. 



The spine term. Since the type space / is finite, we trivially note that (v^,Vp\)^ < oo. It is 
always useful to first focus on the spine term, since we can change the measure with (25) to get 



where, for all ?/ e / 



(32) 



as < — *■ oo and {gtVpX, "^pa)^ = {vx, '^p\) ^ for all t > 0, since rjt is a finite-state irreducible Markov 
chain under with invariant distribution 7r^(y) = v^{y)'^TT{y). It follows that the long term 
the growth or decay of the spine term is determined by the sign of pE\ — Ep\ . 

The sum term. We now assume that pEx — Ep\ > 0. Wc know that under Q\ and conditional 
on knowing 77, the fission times {Sk : fc > 0} on the spine occur as a Poisson process of rate 
(1 + m{r]s))r{ris) with the fc*'^ fission yielding an additional Ak offspring, each A/- being an 
independent copy of A(y) which has the size-biased distribution {{l + k)pk{y) / {l+m{y)) : fc > 0} 
where y = rjSk the type at the time of fission. We also recall from Lemma 2.14 that 

M,{y) := QA(i'(y)) < cx) ^ P{AP{y)) < cx). 

Therefore, if we condition on Qt which knows about (Cs,'7s) at all times < s < t we can 
transform the sum into an integral, use Fubini's theorem and the change of measure used in 
(32): 



fe=i fc=i 



'^l'^ / (l + m(77,))r(7?,)A/,(r,,)«A(r?.)«e«^«=-«^^Ms 



^((1 + m(77,))r(77,)M,(77,)«A(7?.)«e«^«=-«^^^) ds 



v\[y) Jo 

vx (y) pEx - Epx 



where 



p 

hsiy) ■■= QlHfiVs)MM ^{rjs)), f{y) := (1 + m(y))r(y), 
^ Vpx ' 

and kt{y):^^{hv{yYV <t) 

with U an independent exponential of rate {pEx — Epx) > 0. Note that, for all y £ I, hs{y) 
(fMqvl^Vpx) and fct(y) t ^00(2/) as t ^ 00, where (ktVpx^Vpx) = {fMqvl,Vpx)V{U < t) T 
{kooVpx,Vpx) = {rMqV^,Vpx) ■ Then, since Mq{w) < 00 <^=5> P{A{wy) < 00, and / is finite, 
we are guaranteed that fcoo(j/) < 00 for all y G I as long as P{A{wy) < 00 for all w G I. 
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Having dealt with both the spine term and the sum term, we have obtained the upper-bound 

and since Zx{tY is a P-submartingale, we find that 

(pEx - Epx) 

and Zx{t) will be bounded in CP^P^^^y) if we have both pEx - Epx > and P{Ap{w)) < oo for 
all w e I. □ 



3.6.2 Proof of Part 2: 

The earlier proof for BBM goes through with minor modification. Exactly as in the BBM 
case, looking only at the contribution of the spine means that Zx is unbounded in CP{P^'y) if 
pEx — Epx < 0. In addition, letting T be any fission time along the path of the spine. 



Zx{T) > {1 + A{'nT))vx{VT) 



where A{r]T) is the number of additional offspring produced at the time of fission. Then, with 
r>igiy) Q((l + i(y))«) < CX3 ^ PiAP{y)) < oo, 

and so Zx will also be unbounded in CP{P^'y) if mq{y) = oo P{AP{y)) ~ oo for any y G I 

(taking a fission time when also in state y) . □ 



3.6.3 Remarks on signed martingales and Kesten-Stigum type theorems 

In the multi-typed BBM, for each A there will be other (signed) additive martingales correspond- 
ing to the different eigenvectors and eigenvalues obtained from solving (19); the Zx martingale 
simply corresponds to the Perron- Frobenius, or ground-state, eigenvalue Ex and (strictly posi- 
tive) eigenvector vx- Since |u 4- < (|it| -I- Iwl)' < + for all u,w € R, the above proof 
will also adapt to give convergence results for signed martingales. In fact, when there is a com- 
plete orthonormal set of eigenvectors, a Kesten-Stigum like theorem would then swiftly follow 
(for example, see Harris [22] in the context of the continuous-type model of the next section). 

4 A continuous-typed branching-diffusion 

The previous finite-type model was originally inspired by the model that we now turn to, origi- 
nally laid out in Harris and Williams [ ' ] . In this model the type moves on the real line as an 
Orstein-Uhlenbeck process associated with the generator 




with 6 > considered as the temperature, 



which has the standard normal density as its invariant distribution: 

1 12 
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The spatial movement of a particle of type ?/ is a driftless Brownian motion with instantaneous 
variance 

A{y) := ay^, for some fixed a > 0, 



and fission of a particle of type y occurs at a rate 



R{y) := ry^ + p, 



where r, p > are fixed, 



to produce two particles at the same type-space location as the parent (we consider only binary 
splitting). The model has very different behaviour for low temperature values (i.e. low 9), but 
most studies have considered the high temperature regime where 9 > 8r. Also, the parameter A 
must be restricted to an interval (Amin, 0) in order for some of the model's parameters to remain 
in M, where 



Generally, unboundedness in a model's rates is a serious obstacle to classical proofs since they 
often depend on the expectation semigroup of the branching process, and unbounded rates tend 
to lead to unbounded eigenf unctions. Here this is the case, but the existence of a spectral 
theory for their particular expectation operator allowed Harris and Williams to get a sufficiently 
good bound in particular for a non-linear term (see Theorem 5.1 of [20]), and therefore to 
prove £P-convergence of the martingale. Other convergence results for various martingales and 
weighted sums over particles for this model also appear in Harris [22], again using more classical 
methods and requiring 'non-linear' calculations. The spine approach we again adopt here is 
both simple and more generic in nature; requiring no such special 'non-linear' calculations, it 
elegantly produces very good estimates that only involve easy one-particle calculations. 

We use the same notation as previously Xt = {(X„(t), Yu{t)^ : u G Nt} to denote the point 
process of space-type locations in M x M, and suppose that the measures jp^-?' : {x,y) G M^} 
on the natural filtration with a spine {J^t)t>o are such that the initial ancestor starts at {x,y) 
and (Xt, (Ct,77t)) becomes the above-described branching diffusion with a spine. 

4.1 The measure change 

Although there are some significant differences, this model is similar in flavour to our finite-type 
model. There is a strictly-positive martingale Zx defined as 



where v\ and E\ are the eigenvector and eigenvalue associated with the self-adjoint (in £^(7r)) 
operator: 



The eigenfunction vx is normalizable against the /^^(tt) norm, and can be found explicitly as 




ueNt 



Qg + -\^A{y)+R{y). 



vx{y) 




where 




-x/9^ -9{8r + 4aX^), 




The 
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We again define the speed function c\ := —Ex/X, and X{9) < is the unique point (on the 
negative axis) at which c\ hits its minimum - further details are given in Harris and Wilhams 
[20]. We are going to use spines to prove the following result, in which the critical case of A = A 
and the necessary conditions for £''(P)-convergence are new results: 

Theorem 4.1 Suppose that A G (A„iin,0). 

1. Let p G (1,2]. The martingale Z\ is [P) -bounded if both pE\^ E.p\ > and pip^ < V'^Ja- 
In particular, for all A € (A(0),O], Z\ is a uniformly-integrable martingale. 

2. Z\ is unbounded in C^{P) if either pE\ — Ep\ < or pip^ > '(/'jIa- 

3. Almost surely under P , Z\{oo) ~ if X < X{9). 

Once again, for each A < we define a measure Q^'^ on {T,!Foo) via 



dP==^y 

so that with Q\ := Qa|.Foo '^^ have 

Zx{t) Zx{t) 



^T1^2"'«A(r?0e'«'-''^"*, (33) 



dP^y 

The facts arc that under Qa: 



^a(O) vx{v)e^^' 



• the spine diffusion has instantaneous drift arjfX; 



• the type process rjt has generator f (^5- — ^^J/^) ^.nd an invariant probability measure 
tta := (v\,v\J ^^v^Ti, corresponding to a normal distribution, iV(0, 2^)i 

• fission times on the spine occur at the accelerated rate of 2R{rit); 

• all particles not in the spine behave as if under the original measure P. 

4.2 Proof of Theorem 4.1 

Proof of Part 1: Suppose p £ (1,2]. Then using the spine decomposition with Jensen's 
inequality and Proposition 2.13 we find, 

P^^nZ^m < e'^v^iy) (Qr(E «A(^5j'e'^«^"'^^^^") +Qr(«A(?7t)'?e'?^«'-'^^*) 

Assume that pEx — Ep\ > and pip^ < tjj^y^. As seen in Harris and Williams [20], we can do 
many calculations explicitly in this model, largely due to the fact that under Q^^f 

r;, ^ N (e-^^>-'y, - e'^A^p^-) , ^ , . 



2/ipA / V 2^pA, 
and the cigcnfunctions v\ have such simple exponential form. For example, 



pf(^(^/.))-QpA'(e^^'^^"-'^-^''n (34) 
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can easily be seen to be finite and bounded for all s > if and only if — i/jp^ — ~ 

pip'^ — tppx ^ '^i ^^'^ j^^^ ^ readily calculated explicitly. 

In fact, more 'natural' conditions for convergence of the martingales would be that 

{RMqvl,Vpx)^ < oo, {vl,Vpx)^ < oo, and pEx - Epx < 0, 

where Mq{y) := Q(A*(y)) with A the size-biased offspring distribution (here, binary splitting 
means A{y) = 1), and we present arguments below that are more generic in nature, at least in 
terms of adapting to other 'suitably' ergodic type motions and random family sizes. Note, the 
last condition above is related to the natural convexity of E\ and, in our specific model, both 
integrability conditions are guaranteed by ~ ^ 

The spine term. On the algebra Gt the change of measure takes the form 



fx' 



Qt 



exp (^^ RiVs)ds + X{^t - a:) - E't 



which we can use on the spine term to arrive at 

ftix,y) := e^-vx{yWx"[vx{vtre'^^^*-'''=^') = e^'-^Vpxiy) gt{y)e-^^'''-''-'^' (35) 

with gt{y) '■= ^px {^^{'Ht) / Vpxivt)^ ■ Under the assumption that pip^ < tppx, it easy to check 

that (v^x^Vpx)^ < oo, that is v^^/vpx G L^{TTpx) from which it follows that gt G L^{TTpx) for all 

t > Q. Since 77 has equilibrium Hpx under QpAi we find {gtVpx,Vpx')^ = {vxjVpx}^ < 00 and 

gtiy) {vl,Vpx)^{vp\,Vp\)^^ < 00 as t ^ 00 for all y e M. 

We also note that since gt G L^{'Kpx), we have ft G L^liTpx) where tt^ := (1,171)^ v^i: and 
then 

/ ^px{y)h{x,y)dy ^ e'"^i<l^e-^pE.-E^^)t_ 

The sum term. Note that under the parameter assumptions we have (^Rv^,Vpx)^ < 00. As 
for the finite-type model the fission times Su on the spine occur as a Cox process and therefore 

gtix,y) := e'^vxiy)Qr [Y. «A(^5j«e«^««" -'^^^") 

h,x{y) ^*Q^'^(2i?(7y,)«A(?7.)^e«'«=-'^^^) ds 



e'^'vxiy) / (2i?(7y,)^(77,))e-('^^^-^-^)Ms 
Jo ^ Vp\ ' 

eP^''vpx{y)h{y) 



where 



kt{y):= f /j,(y)e-(^'^^-^-^)Ms, /i,(y) := Q;;,^ (2i?(r;,) ^(77,)) 

Jo ^ "f^pA ^ 



and ht, kt G L^{Trpx)- Note, kt{y) t ^"00(2/) G L^{Trpx) as i ^ 00 where 
{ktvpx,vpx), (l-e-(P^--^^-)^) ^ {2Rvl,Vpx), 

— ^ = {2Rv'^,Vpx}^ ^ , , T -— -=r-T {kooVpx,Vpx)^ < 00. 

{vpx,Vpx)_ A F {pEx-Epx) [pEx-Epx) ^ fi-w 
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Note, kt € L^{Trp\) implies gt S L^{Trp\), with an explicit calculation again possible. 
Bringing together the results for the sum and spine terms, we have an upper bound 



and hence the submartingale property reveals that 

for alH > and all y eR. □ 



Proof of Part 2: As there is family sizes pose no problems, we need only dominate the 
martingale by the spine at fixed time t: 



v\(y) P \Vpx ) 

and Z\ is therefore unbounded in D'{P'^) if either -pEx — Ep\ < or (u^, Vp\) ^ = oo. □ 



Proof of Part 3: The proof that we have seen in the finite- type model will work here with 
little change: under Q;,, the spatial motion is 



and the type process 77^ has invariant distribution A^(0, 5^)' whence t ^£_t ^ Aa0//iA = i?^. 
Therefore it follows that under the diffusion + c^t drifts off to — 00 if A < A(0). When 
A = A, it is also simple to check that + c\t is recurrent, so has liminfj^j + cxt} — —00. 
Whence, in either case, bounding Zx below by the spine's contribution as done before, we have 
Zx{t) > vx{r]t)e^'^^*'^'^^^^ and since vx > and r]t recurrent, we see that limsupj^o^ ^A(i) = 00 
fx' 



almost surely under Q^'^. □ 
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